Applications of Temperley-Lieb algebras to Lorentz lattice gases 
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Abstract 

Motived by the study of motion in a random environment we introduce and investigate 
a variant of the Temperley-Lieb algebra. This algebra is very rich, providing us three 
classes of solutions of the Yang-Baxter equation. This allows us to establish a theoretical 
framework to study the diffusive behaviour of a Lorentz Lattice gas. Exact results for 
the geometrical scaling behaviour of closed paths are also presented. 
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The Temperley-Lieb(TL) algebra arosed in the context of statistical mechanics and it was 
used to map the g-state Potts models into the 6- vertex model M. The relevance of this algebra 
to the theory of two-dimensional solvable models was further elaborated by Baxter || and the 
first detailed study of its mathematical structure was given by Jones ||. In recent years, the 
interest in TL algebras is widespread in many branches of physics and mathematics, including 
quantum spin chains f|, |], conformal field theory |J and knot theory J^, |8| just to name a few. 

In this letter, motived by physical considerations, we introduce and investigate a variant 
of the TL algebra. Our physical motivation is to describe the diffusion behaviour of a Lorentz 
lattice gas, consisting of a dilute gas of particles moving in a random array of scatterers || 
. In the model we shall consider here, a particle moves along the bonds of the square lattice 
and its trajectory is changed when it hits scatterers which are either (right and left) mirrors 
or rotators randomly placed at the sites of the lattice [0, [TIJ. This is perhaps the simplest 



nontrivial "scattering rules" one can think of and they are illustrated in Fig.l . To describe 
the statistical behaviour of these paths in a more general way, we assign a fugacity q to every 
closed trajectory. This then defines a loop model and its partition function is given by 

Z= wl a wl b w n c c w n d d q* paths (1) 

scatter configurations 

where w a (w c ) and w (wd) are the Boltzmann weights probabilities of right and left mir- 
rors(rotators) and rij(i = a,b,c,d) are the the number of weights in a given configuration. 
Setting q = 1 we recover the Lorentz gas model. 

In the absence of rotators (w c — Wd — 0), the partition function (1) can be interpreted as 
a graphical representation of the g 2 -state Potts model, and therefore the algebraic structure 
underlying the mirror collision rules is the standard TL algebra. This is no longer valid if 
rotators are present, because a particle hitting a rotator does not behave exactly the same 
when it hits a mirror. In fact, the algebra that mimics the behaviour of rotators is as follows. 
We associate to right and left rotators the elements Ri and Li, acting on sites i and i + 1 
of a quantum spin chain of size L. We find that Ri and Lj are generators of the following 
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associative algebra 

R* = qRi, L? = qLi, [Ri,Li} = (2) 

LiRi±\Li = Li, RiLi±iRi = R iy i? i±1 ] = [Lj, Lj ±1 ] = (3) 

[L u Lj] = [R u Rj] = [Ri, Lj) = \i - j | > 2 (4) 

These relations can be seen as two commuting TL algebras whose generators alternate 
between R t and Li depending on the parity of the sites. This algebra has to be read in 
conjunction with the TL operator Ei = RiLi and the identity Ii, which in the Lorentz gas 
corresponds to left and right mirrors, respectively. It is interesting to note that this algebra 
also allows us to define the braid operator 6, = Ri + Li — e 9 Ei — e~ e Ii and its inverse b^ 1 = 
Ri + Li — e~ 9 Ei — e d Ii, where the parameter 9 is related to the fugacity by q = 2 cosh(6 l ). The 
real surprise is the fact that the braid-monoid operators &j, 6" 1 and Ei satisfy the Birman- 



Wenzel- Murakami (BWM) algebra [11 1. More precisely, this connection occurs when the two 
independent parameters y/Q and c of the BWM algebra (we are using notations of ref.||) lie 
on the curve \fQ = 4 cosh 2 (9) and c = — e 38 . However, in the context of the Yang-Baxter 
equation, we find that the "rotator" algebra (2-4) is richer than the BWM algebra(on the 
above curve), since the former gives us an extra integrable manifold. 

To make further progress we look for possible integrable manifolds for the Boltzmann 
weights w a , Wb, w c and w,i- This not only allows us to establish a theoretical framework to 
study the Lorentz lattice gas but also prompts us to make exact predictions for its diffusive 
behaviour. The basic idea is to search for solutions of the Yang-Baxter equation by Baxterizing 
the Ansatz w a Ii+WbEi+w c Ri+WdLi with the help of the algebraic relations (2-4). Here we omit 
technical details, presenting only the final results. We find three different classes of integrable 
manifolds and they are given by 

(/) w a (X) = SmH6 ~ X \ w b (X) = l/w a (X), w c (X) = w d (X) = 1 (5) 
smn(A) 

. i sinh(# — A) cosh(6 l — A) 
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2A 

(III) w a (X) = w b (X) = q2 ^l_ e2X , Wc(A) = 1, w d (X) = (7) 

where A is the spectral parameter of the Yang-Baxter equation. The third solution also admits 
w c (X) = and u^(A) = 1, since the algebra (2-4) is invariant under the exchange of left and 
right generators. 

In order to interpret these solutions it is more appealing to work with a specific representa- 
tion for the generators Ri and Li. In the language of quantum spin chains, we find that these 
generators can be written as 

ft = o t r- +i + ."ri, + ^(1 - „>r> +1 ) + ^fl(r Ul - <) (8) 

U = r?«- +1 + n.i, + SS*W (1 - r^ ?+1 ) + - r< ) (9) 

where {erf, of} and {r,- , rf } are two commuting sets of Pauli matrices acting on the sites of a 
lattice of size L. 

Considering this representation, it is not difficult to see that solution I behaves similarly as 
two decoupled 6-vertex models and therefore the underlying symmetry is based on the D 2 Lie 
algebra ||12|| . The second solution corresponds to a fine tuned coupled 6-vertex models through 
their (total) energy-energy interaction and its quantum group symmetry, after a canonical 
transformation, can be related to the twisted A\ Lie algebra ||12|| . The third solution is clearly 
asymmetric in the spin variables and therefore cannot be interpreted in terms of the BWM 
algebra. This solution has also a number of unusual properties. For instance, the asymptotic 
braid limits _R(±oo) are not invertible and for the special values q = ±1 we find a connection 
to a singular point (t = in the notation of |ll[) of the Hecke algebra. In the context of 
spin models, we have verified that the Boltzmann weights (7) can be related to an integrable 
non-self-dual manifold of generalized Potts models introduced by Domany and Riedel |n 
to model the adsorption of molecules on a crystal surface. 

We now turn our attention to the Lorentz lattice gas. For q = 1 only the first two solutions 
turns out to be physically meaningful, since all the Boltzmann weights can be made positively 
definite. Thus, we shall concentrate our efforts to study the critical behaviour of these two 
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solutions with the expectation of making predictions for the geometrical scaling behaviour of 
closed particle paths in the Lorentz lattice gas. To this end, we study these models with 
generalized boundary conditions, in order to assure that closed loops on the cylinder pick up 
the correct Boltzmann weights. In general, these are twisted boundary conditions defined by 

= rt +1 = e^r* al +1 = of, r[ +1 = r{ (10) 

where ip\ and ip2 are arbitrary angles. 

The critical behaviour of solution I is that of two decoupled XX Z spin chains and conse- 
quently the results for the central charge and anomalous dimensions can be read off directly 
from previous work in the literature ||. This system is critical in the regime q G [2, —2] and it 
is more convenient to parametrize the fugacity by q = 2cos(7). It turns out that the effective 
central charge behaviour is 



while the conformal dimensions are given by 
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-y-mi,m2 \ " 



i=l 



4x p (7) 



(12) 



where n» and m« are integers representing the spin-wave and the vortex excitations of two 
decoupled Coulomb gas having the same radius amplitude 2^(7) = '^ L . The finite-size correc- 
tions (11,12) are measured relative to the ground state of two periodic XX Z models 

The second solution corresponds to two nontrivialy coupled XX Z models and it is possible 
to show that its Bethe Ansatz solution is formally related to that of the A$ vertex model [ITSl . 



The critical properties of the latter model has only been partially studied in the literature |16 
and in a region that exclude the Lorentz lattice gas itself. We remark that here we have to 
perform the calculations in the presence of the seam ip\ and ip2, since they plays a crucial role in 
the underlying critical behaviour. We find that the eigenvalues of the Hamiltonian associated 
to the second solution are given by 

En(L) = ± 2s -° 2 w (13) 

j=i cos(7) — cosh(A^ ) 
4 



and the corresponding Bethe Ansatz equations are 
sinh(A5 1) /2 - i7/2) 



^ sinh(Af } /2 - AgVg - ij) u siDh(A^ - Aj 2) + ij) ._ 



tA sinh(Aj 1) /2 - X^/2 + 27) i 1 ! smh(Af - A< 2) - z 7 ) 



sinh(A5 1) / 2 + *7/2) 

(14) 



L\ sinh(Af } - Af - *2 7 ) ti sinh(Af - A« - z 7 ) 
The existence of the Bethe Ansatz solution allows us to calculate, besides the thermody- 
namic limit, the dominant finite-size corrections for the eigenvalues of the Hamiltonian. For 
a conformally invariant system, these finite-size effects can be directly related to the central 
charge and scaling dimensions |T7|, providing us a way to study the universality class of 



solution II. An essential step here is to determine the nature of the Bethe Anstaz roots gov- 
erning the ground state properties. In the regime we are interested, i.e. near 7 ~ 7r/3, we find 
a mixture between one-string and two-string type of solutions. More precisely, the complex 
roots structure are given by 

\f = ± *(vr/2 - 7) + 0(e- L ), Xf = £f + «r/2 (16) 

where t^f* are real numbers. 

In order to determine the finite-size corrections for the ground state, we have numerically 
solved the Bethe Ansatz equations for several values of the lattice size up to L ~ 40. In table 
1 we present our estimates for the effective central charge for general boundary conditions. 
Surprisingly, the behaviour is precisely the same of that given in formula (11) which suggests 
that the coupling between the two XX Z models becomes asymptotically irrelevant near q = 
1. To give extra support to this scenario we have also investigated the low- lying spin- wave 
excitations. In table 2 we show some estimates for the exponents X®'q and Xi'i and they 
are in accordance with the conformal dimensions (12). More generally, we verified that the 
critical behaviour for role region q G [0, \/2) is given in terms of two decoupled Coulomb- 
gas models. However, in the regime q G [v2, 2] the coupling between the two XX Z models 
becomes relevant, and the criticality is governed by a c = 3 conformal field theory. We note 
that in the first branch all the Boltzmann weights probabilities can be positively defined while 



in the second one either Wf, or w a is necessarily negative. This is perhaps the physical reason 
behind these two different anti ferromagnetic critical behaviour. Finally, we remark that our 
results complement and correct an early calculation in the literature (16| . 

We now have the basic ingredients to study the scaling properties of closed trajectories 
in the Lorentz lattice gas model. From our previous analysis we conclude that the critical 
behaviour, for both manifolds I and II, is given in terms of two decoupled Coulomb-gas in the 
presence of the background charges ipi = 4>2 = 27 = 2ir/3. Here we are interested in correlators 
that measure the probabilities that two points on the lattice separated by r belong to the same 
loop [T9|]. For large distances we expect the algebraic decay 



< $;(r)$i(0) >~ r 



-2X, 



(17) 



where X\ is the scaling dimension of the conformal operator $; (r) . Equipped with equations 
(11) and (12) we can now calculate the dimensions X\ for the q = 1 model. The probability 
that / = ni + n 2 loop segments meet in the neighbourhood of two points in the lattice is 
associated to the conformal dimensions of the spin-wave excitations n\ and n 2 with null vortex 
charges m x 



m 2 = [|T^, 20 1. Due to the background charges, we also have to subtract the 
variation of the true ground state for the Lorentz lattice gas. We conclude that the value of 
these dimensions is 



X., 



ni,n2 



E 

1=1 

2 

E 

i=i 



71? (1 



(1/3) 2 | 
4x p (tt/3) 

<5n;,o) 



5, 



rii.OJ 
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The scaling behaviour of single paths is governed by the lowest conformal dimension X 10 = 
X 0i i = 1/4, predicting a fractal dimension []2l] df = 2 — 1/4 = 7/4. This is the same value 



found in numerical simulations of a fully occupied lattice of either mirrors or rotators [21, 22 



Therefore, our result brings an extra theoretical support to the fact that the scaling behaviour 
of single paths does not depend whether the scatterers are mirrors or rotators. However, the 
situation changes drastically when we consider the probability for multi loops I > 2. In fact, 
in the case of mixed rotator-mirror model the second lowest dimension occurs in the sector 
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rii = n% = 1, and the corresponding fractal dimension is df = 2 — 1/2 = 3/2. This value is 
double of that expected df = 2 — 5/4 = 3/4 when we have only mirrors |20|| . In general, formula 
(18) leads us to conclude that the multi loops scaling behaviour is sensible to the details of 
the scattering mechanism. Finally, we remark that our findings are valid for rather distinct 
manifolds, suggesting that w c = Wd is indeed a critical surface in the Lorentz gas model. 

In summary, we have introduced a rather rich variant of the Temperley-Lieb algebra which 
generates three classes of solutions of the Yang-Baxter equation. We have applied the first 
two of them to study the diffusion behaviour of a Lorentz lattice gas whose scatterers are 
either mirrors or rotators. The third manifold seems to be relevant to the study of phase 



transition in adsorbed films [ 14 1, and we hope to investigate its critical properties in a future 



publication. Finally, we remark that after this paper has been written we noted the very recent 



work P3J in which the third manifold and its generalizations were obtained via Baxterization 



of Fuss-Catalan algebras. 
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J 


<Pl — — TV/ O 


l/l-i 1/ln. TV / A 

ipi — lp2 — 71/ 41 


l/l-i TT / ~\ 1/li-. TT / /I 

ty^l — 7T/0, ^2 — 71/ 4 


o 




1.55667 


1. 78890 


1.67253 


lb 


l.oloo I 


l.boUl ( 


l.boUl < 


24 


1.51021 


1.73159 


1.62083 


32 


1.50683 


1.72737 


1.61705 


40 


1.50506 


1.72516 


1.61507 


Extr. 


1.5001(2) 


1.7186(2) 


1.6092(2) 



L 


*S(7 = 7T/3) 


*J°( 7 = tt/3) 


X 1 ; ( 7 = tt/3.5) 


X 1 ; 1 ( 7 = tt/2.5) 


8 


0.344610 


0.640435 


0.404331 


0.599395 


16 


0.337923 


0.653822 


0.392503 


0.599587 


24 


0.336174 


0.658341 


0.387451 


0.599756 


32 


0.335383 


0.660548 


0.384344 


0.599839 


40 


0.334934 


0.661844 


0.382155 


0.599890 


Extr. 


0.33330(2) 


0.66665(2) 


0.3575(3) 


0.6000(2) 



Table 1 : Finite size sequences for the extrapolation of the effective central charge for 7 = 
tt/3. 

Table 2 : Finite size sequences for the extrapolation of the spin-wave anomalous dimensions 
X°o and X®'® for periodic boundary conditions. 
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